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A warped m = 2 water maser disc in V778 Cyg? 
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ABSTRACT 

The silicate carbon star V778 Cyg is a source of 22 GHz water maser emission which was 
recently resolved by MERLIN. Observations revealed an elongated iS-like structure along 
which the velocities of the maser features show a linear dependence on the impact parameter. 
This is consistent with a doubly-warped m = 2 disc observed edge-on. Water masers and 
silicate dust emission (detected by IRAS and ISO) have a common origin in O-rich material 
and are likely to be co-located in the disc. We propose a detailed self-consistent model of 
a masing gas-dust disc around a companion to the carbon star in a binary system, which 
allows us to estimate the companion mass of 1.7 + 0.1 M , the disc radius of 40 ± 3 AU and 
the distance between companions of ~ 80 AU. Using a dust-gas coupling model for water 
masing, we calculate the maser power self-consistently, accounting for both the gas and the 
dust energy balances. Comparing the simulation results with the observational data, we deduce 
the main physical parameters of the masing disc, such as the gas and dust temperatures and 
their densities. We also present an analysis of the stability of the disc. 

Key words: circumstellar matter - masers - stars: AGB and post-AGB - stars: carbon - stars: 
chemically peculiar - stars: binaries 



1 INTRODUCTION 

Carbon stars with amorphous silicate dust features at 10 and 18 
Hm (silicate carbon stars) which typify O-rich dust envelopes 
lLittle-Marenirtfl986l) . were first discovered with the Infrared As- 
tronomical satellite (IRAS). A distinguishing feature of this class 
of carbon stars is a mixed chemistry, i.e. simultaneous presence of 
O- and C-rich material around them. The Infrared Space Observa- 
tory (ISO) revealed that this mixed chemistry phenomenon is not 
unique . About 20 silicate ca rbon stars are known in our Galaxy 
IChen et alll999tlSzczerbal2002h . 

The best known example of this class is V778 Cyg. Infrared 
observations show that the shape and intensity of the silicate fea- 
tures in the spectr um of V77 8 Cyg did not change during a period 
of about 14 years lYamamura et all2000l hereafter Y00). This in- 
dicates that the O-rich material is located in some stable configu- 
ration. Y00 proposed that it is stored in a disc around an invisible 
(most probably main-sequence) star, whi ch is the companion of the 
carbon star in a binary syste m (see also M orris et alJ H 1987 : Evans 
1199(1 fearnbaum et alJll99ll) . V778 Cyg also shows H 2 and OH 
maser emission. Single dish monitoring of water maser emission 
during 12 years shows a very small velocity shift (< 0.5 km s -1 ) 



of the main components in the maser spectrum lEngels & Leinertl 
I1994 1 . 

In October 2001 water maser emissio n from V778 Cyg wa s 
imaged using five telescopes of MERLIN ISzczerba et alJ l2006). 
The position of the maser complex was found to be displaced by 
~ 190 mas from the position of the C-star measured 10 years previ- 
ously, which is unlikely to be a result of the proper motion. Instead, 
this can be considered as an additional evidence that the O-rich ma- 
terial is stored around the secondary ISzczerba et al.l2006h . 

The MERLIN data show that the masers are distributed in an 
<S-shaped figure. In this paper, we construct a model of the spatial 
and velocity distribution of the water masers around V778 Cyg. We 
develop a self-consistent description of the physical mechanism of 
the water maser pumping together with an explanation for the ob- 
served infrared silicate dust emission. We also determine the main 
parameters of the system. 

We present the main observational data in Section |2| We 
model the 22-GHz water maser spatial and velocity distributions 
and the infrared dust emission as well as the physical mechanism 
of maser pumping in Section|3| We discuss our results in SectionfH 
and conclude in Section|5| 
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Figure 1. Water maser distribution observed 
towards V778 Cyg ISzczerba et"ai]l200fil) . 
The coordinates are measured relative to 
the position of the strongest feature at -15 
km s _I at (a,S) = (20 h 36 m 07 s .38327, 
60°05'26'.'024). The solid curve shows the 
best fit to a model of a doubly-warped ring 
of radius 40 AU (at the assumed distance of 
1.8 kpc) which is thin in the radial direction; 
other parameters are presented in Table 1. 
The ring centre is at the origin of the coor- 
dinate system. The insert shows the MER- 
LIN spectrum of the H2O maser emission. 
The dotted curve shows the spectrum magni- 
fied by a factor of four to enlarge the weak 
features. 
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2 DATA 

The silicate carbon star V778 Cyg is located at a distance of about 
d = 1.8 kpc, whi le the carbon star luminosity is L, m 5.6 10 3 L Q 
jBenson & Little-Mareninll996l) . Observations of V778 Cyg at 22 
GHz by MERLIN revealed that the water maser f eatures form an 
<S-sha ped structure with an angular size of 18 mas (Szczerba et al. 
2006; see Fig.0. The observed line-of-sight velocities of masing 
blobs depend almost linearly on the impact parameter measured 
along the major axis of the structure as shown in Fig. [2] The total 
spectrum of water maser emission (see insert in Fig.0 consists of 
four major components near -22, -19, -17 and -15 km s -1 . The 
last component is the strongest, with a central peak of about 15 Jy, 
which corresponds to a lower limit to the brightness temperature of 
T b = 6 10 s K ISzczerba et all2006h . 



3 MODELLING 



ISzczerba et alj feOOot) show that the structure and velocity field 
of the maser features (presented in Figs. Q and are consistent 
with a Keplerian disc, observed almost edge-on, and present a sim- 
ple quantitative interpretation of the observational data. The three 
maser spectral components at -19, -17, and -15 km s _I (see 
Fig. have been detected in every observation during the whole 
period from 1988 to 2002. They are symmetric around -17 km s 
in velocity space. They are also almost equally displaced on the 
sky. By analogy with water masers in AGNs, this has prompted the 
interpretation that these are the systemic and the high velocity com- 
ponents that appear towards the centre and in the wings o f a Kep- 
lerian disc (see e.g. lMivoshi et alll995tlBragg et all200 0). In that 
case the velocities of the individual features making up the systemic 
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Figure 2. Dependence of the velocities of the maser features on the impact 
parameter measured along the major axis of the observed structure, relative 
to the brightest feature at -15 km s - ', which is assumed to be the systemic 
velocity V sys . The solid curve represents the model velocity dependence 
corresponding to a Keplerian disc with the inner and outer radii of R m = 10 
and R out = 40 AU, with V(b) computed from dr m /db = 0, where r m is 
given by Eq. The linear dependence V(b) = VK(R ut)b I Rout is shown by 
a dotted line for comparison. The pure Keplerian rotation law V(b) oc 1 / yfb 
is shown by a dashed curve. 




Figure 3. The optical depth in a masing line r ra (in units ArR ollt ) forf = Vn(R ut) I AVd = 11.7 and the ratio of the inner-to-outer disc radii R m IR oa t = 1/4. 
The parameters are chosen to represent a disc around a 1.7M Q star, with the outer radius R out = 40 AU and the gas temperature of 300 K. (a) Dependence 
on the impact parameter b for the velocities V = V/AVd = (dotted curve), 5.7 (dashed curve), 11.4 (dot-dashed curve), 13.3 (triple dot-dashed curve). 
(b) Dependence on the velocity of maser features for h/R out = (dotted curve), 0.3 (dashed curve), 0.6 (dot-dashed curve), 0.9 (triple dot-dashed curve). 
Solid curves show the maximum r m as a function of b and V (see Section lTTl . 



component should show a linear dependence on the impact param- 
eter V(b) oc b jGrinin & Grigor' evil 19831: IWatson & WyidM200Cl) 
which indeed is the case (see Fig.[2}- This linear dependence could 
be interpreted as solid-body rotation. However, such an explana- 
tion is highly unlikely, because this requires a self-gravitating disc 
having a mass comparable to or even larger than the central mass. 
In such a case, the number density in the disc would be too high 
for the maser to operate as the population levels of water molecules 
would become thermalised by collisions with hydrogen. 

On the other hand, the Keplerian disc interpretation implies 
a velocity at the outer radius of Vk ~ 2 km s and a disc radius 
of if - 11 AU. The component at -22 km s _1 , which appeared 
in the spectrum in 1992, causes a problem for this model. It lies 
farthest from the disc centre but has too large a velocity for a simple 
Keplerian disc. 

The <S-shape is l ikely to result from the curving of the disc 
ISzczerba et all TToOot) . However, the observed displacement of 
maser features at -19 km s -1 relative to those at -17 ± 1 km s -1 
is very large and is comparable to the estimated disc radius, which 
would be inconsistent with a simple model of a slightly warped 
masing dis c, widely accepted, for example, for NGC 4258 (e.g. 
iMivoshi et all 1 99a) . The disc radius thus could be larger than the 
assumed 1 1 AU (corresponding to 6 mas), and the maser emission 
at -22 km s -1 could still come from the disc. 

The stability of t he dust spectrum (Y00 ) and of the water 
maser spectral features l En gels & Leinertfl9 94) prompts us to sug- 
gest that they have a common origin in the O-rich material in the 
disc. This disc is unlikely to be located around the carbon star, be- 
cause the dust would be blown away by the stellar radiation at the 
distance corresponding to the observed disc size. It is thus natu- 
ral to accept the hypothesis of Y00 and consider a disc around the 
companion of the carbon star in a binary system. 



3.1 Masers from a Keplerian disc 

Let us first consider a flat thin disc in Keplerian motion. 
The theory of masers from Keplerian discs is w ell developed 
iGrinin & Grigor' evl Il983l : IWatson & Wyldl l200fj) . Masers are 
bright when the line of sight crosses the disc at a grazing angle, 
close to edge-on. Due to the velocity gradient, the maximum maser 
coherence length is reached towards the disc centre (where the mat- 
ter moves perpendicular to the line of sight) and near the outer parts 
(where the gas motion is almost parallel to the line of sight). The 
strongest spectral features thus appear at the systemic velocity and 
at velocities shifted towards the red and blue by almost the Keple- 
rian velocity near the outer edge of the masing disc. 

The optical depth in a masing line, integrated along the line 
of sight s, as a function of the impact para meter b and veloc- 
ity V is given by (s ee, for example, IGrinin & Grigor' evlll983t 
IWatson&Wvlj2 000) 



r m (b, V) = 




where k is the maser opacity, = ^kW/^is the Keplerian an- 

gular velocity at radius r = Vs 2 + b 2 , AVd = (2kT /mH,o)'' 2 is the 
thermal velocity for water molecules, and T is the gas temperature. 

The upper integration limit -Jr^ — b 2 is determined by the outer 

disc radius R out . The lower limit is either zero or -^ 2 n - b 1 (corre- 
sponding to the inner disc radius Rj n ). Equation Q gives half of the 
total optical depth along the line of sight. The optical depth can be 
expressed as a function of dimensionless variables b' = b/R out and 
V = V/AV D : 
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Table 1. Model parameters 



Rout" ' b <P c & d A" V sy / M s g 

AU deg deg deg km s _I M Q 

40 ± 3 84.5 ±1 -8 ±2 4±2 0.23 ± 0.01 -15 1.7 ±0.1 



the masers are produced in a thin curved disc of a shape described 
by the following expression 

z(R, <p) = A(R) R sin[m(0 - ©)], (3) 

where R is the cylindrical radius, (f> is the azimuthal angle and z 
is the deviation from the plane XY, where the Z-axis is chosen in 
the direction of the disc angular velocity. A is the amplitude of the 
shear wave (which we assume to be constant), and m and O are 
its azimuthal number and phase, respectively. The projection of the 
disc on the sky depends on the inclination angle ( (the angle be- 
tween the line of sight and the Z-axis) and on the position angle of 
the line of nodes 0. 

Theoretically, it is rather easy to produce a warp with az- 
imuthal number m = 1 in a disc in a binary system, using tidal 
effects. The theory of such discs is quite well developed (see, for 
example, Papaloi zou & Terauemll 19951 and references therein). In 
that case, rings at different radii can have different inclination an- 
gles and the observed spatial distribution in Fig.Qcould be repro- 
duced by two such rings: one ring of a larger radius reproduces 
maser features at velocities V > -15.9 km s _1 and V < -18.5 
km s~', while another ring of a smaller radius produces features at 
-18.3 km s -1 < V < -16 km s" 1 . However, the observed linear 
increase of maser velocities with the impact parameter means that 
all maser features lie on very close orbits. This makes the m = 1 
warp interpretation improbable. 

On the other hand, we find that the observed <S-like structure 
can be well described by a single orbit with azimuthal mode m = 2. 
We discuss possible reasons for origin of such warp in Section l4"2l 
We fit the spatial distribution of maser features and obtain the set of 
parameters which are given in Table 1 (parameters R = i? out , i, 0, 
and A). For illustration we present a three-dimensional image of 
the disc taking R out = R = 40 AU, assuming that A is constant and 
Rir, = 10 AU, in Fig.|4ja)- At large impact parameters the coherent 
length for maser amplification is large in case of the flat disc. For 
the warped disc shown in Figure|4| the line of sight passes through a 
large portion of the disc at large negative values of X, and the maser 
is bright. However, at large positive X (see the extreme dotted curve 
in Fig.[4j> corresponding to X = 15 mas) it crosses the disc plane 
at a steep angle reducing thus the amplification length. At small 
negative X (see the thick solid curve 2), the line of sight passes 
through the whole disc, while at small positive X (dotted curves) it 
passes though one part of the disc only. This would explain why we 
only detect masers from one side of the disc. We explain this more 
fully in the next Section. 

3.3 Maser velocity distribution 

In a doubly-warped disc with m = 2, the velocities are not precisely 
described by the Keplerian rotation law. However, because the warp 
is not very large, we neglect possible deviations from Keplerian 
motion. To describe the velocity distribution of maser features, we 
fit the central mass M s , while keeping all other parameters fixed at 
the values deduced from the geometry in Sect. 13. 21 

In Section l3~2l we showed that the observed dependence of 
the velocity on the impact parameter can be well described by a 
theoretical V(b) dependence for a Keplerian disc with M s = 1.7 M Q 
and other parameters from Table 1 (see solid curve in Fig.[5). In our 
interpretation the strongest maser feature at -15 km s _I coincides 
with the systemic velocity (gas moving perpendicularly to the line 
of sight) and is projected towards the disc centre b = 0. The features 
at « -17 km s (i.e. -2 km s on Fig.|5J corresponding to the 
impact parameter b ^ 0.2/? out show a linear dependence close to 



" The outer disc radius. * The inclination angle of the disc plane. c The 
phase angle (see eq. 0). d The position angle of the line of nodes. e The 
amplitude of the shear wave (see eq. ' Systemic velocity. s The mass 
of the companion. 

where f and s' are the corresponding distances measured in units 
of i? out , and £ = Vk^ouiVAVd is the ratio of the Keplerian velocity 
at the outer disc edge to the thermal velocity. At constant opacity 
k and temperature T, the integral in Eq. J2j depends only on £ and 
on the disc geometry expressed through the ratio of the inner to the 
outer disc radii R m /R OU f 

Fig. shows the dependence of the optical depth in a mas- 
ing line r m on the impact parameter b and on the maser velocity 
V (the figure is identical if we substitute v' — > —V',b' — > —b'). 
Panel (b) shows that r m as a function of velocity has two peaks: 
near zero (systemic) velocity, where the impact parameter b ~ 0, 
and near the Keplerian velocity at the outer disc radius V ~ f . In 
the latter case the dominant rays emerge at an impact parameter 
b x 0.9R out (for f = 1 1.7, which is considered here for application 
to V778 Cyg). In general, for a fixed b, the velocity spectrum has 
two maxima: one corresponds to almost solid-body rotation with 
V(b) as VK(R ut)b/R 0U t, another corresponds to the Keplerian rota- 
tion with V(b) =5 V K (fe). The first one is much stronger at small 
impact parameters, while the latter dominates at larger b. 

For a given velocity V, the optical depth as a function of the 
impact parameter b has a single maximum, as in Fig.|3Ja). If emis- 
sion at any given velocity i s only detected at a single spatial position 
(e.g. lSzczerba et all2006t) . V(b) can be evaluated by fixing V and 
finding the value of b for which r m reaches the maximum, i.e. by 
solving the equation dr m /db = 0. This gives the relationship shown 
by the solid curve in Fig.|5| It has a solid-body rotation part at small 
b and a Keplerian part at large b. We also note that V(b) is almost 
independent of £ and R in /R out . 



3.2 Disc geometry 

The longest maser amplification length is reached when radiation 
passes through the main body of the disc or through its outer parts. 
The strongest spectral line for Keplerian discs (with small ratio of 
inner-to-outer radii) observed edge-on appears at the systemic ve- 
locity and its spatial position coincides with the projection of the 
disc centre. We therefore associate the strongest component in the 
observed maser spectrum of V778 Cyg (chosen as the origin of co- 
ordinates in Fig.0 with the systemic component coinciding with 
the projection of the disc centre and fix the systemic velocity V sys at 
-15 km s~', which is consistent with a value of V svs ob tained from 
the ob servations of OH masers (-15.7 km s '. fLittle-Marenin et all 
Il988t -15.1 km s '. lBarnbaum et alll99ll) . We also expect to see 
strong maser features near the projection of the outer radius of the 
masing disc at velocities shifted from V sys by the Keplerian velocity 
at this radius. 

We would like to describe the positions of all maser features 
(Fig.0 and their velocities (Fig.|2l simultaneously. We assume that 
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Figure 4. (a) Three-dimensional image of the mid-plane of the masing disc with m = 2 and of the maser spots in V778 Cyg. The line of sight lies in the YZ 
plane and makes an angle of 5.5° with the F-axis. The coordinates of the disc centre are (0,0,0). The grey (green) spots show positions of masers. 
(b) Cross-sections of the disc mid-plane for X = - 13, -8, -5, 0, 2, 4, 15 mas. Solid and dotted curves correspond to negative and positive values of X. The 
strongest masers at -22, -17 and -15 km s appearing at X a -13, -5 and are shown as spots in panel (a) and marked by 1, 2 and 3, respectively. The disc 
cross-sections are these X are shown by thick solid curves in panel (b). The dashed line with an arrow shows the line of sight. 



V(b) as Vn(R aui )b/R out (shown by the dotted line). If the disc were 
flat, these features would be weaker than the systemic feature at 
b = 0, because the optical depth at this position is about 40 per 
cent of that at b = for the assumed gas temperature T = 300 K 
(see Fig. [3}. However, because of the warping, the line of sight at 
b a 0.2R aui passes almost through the whole disc (see curve 2 on 
Fig.|Ub)), whereas the amplification at b = takes place through 
only half of the disc (again because of the warp). This results in a 
strong increase of the maser power at these intermediate velocities 
and impact parameters. 

The weakest detected component at -22 km s -1 (i.e. -7 km s _I 
in Fig.|2j lies exactly at the velocity-impact parameter relationship 
computed for a Keplerian disc, which at large b follows a simple 
Keplerian rotation law V(b) oc 1/ sfb (the dashed curve in Fig.0. 
The position and velocity of this component are consistent with the 
expected theoretical peak of the maser emission at b « 0.9R out and 
V * V K (R out ) (see Fig.0. 

We now note that for a flat, symmetric disc both blue-shifted 
and red-shifted components should be observed around the sys- 
temic velocity. However, only blue-shifted features are visible 
(Fig.Q. In a warped disc the line of sight cuts the near and far 
sides of the disc at different angles. The cross-sections of the cen- 
tral disc plane are shown in Fig.EJb). The solid and dotted curves 
correspond to the projections where the maser amplification would 
give the blue- and the red-shifted components, respectively. It is 
clear from the discussion in the previous section that the ampli- 
fication lengths for the blue-shifted maser components are larger, 
which would explain the absence of the red-shifted components in 
the maser spectrum. 



3.4 Dust emission 

The infrared spectrum of V778 Cyg obtained by 750/SWS and 
IRAS/LRS shows strong silicate features at 10 and 18yum (see Fig.|5| 
and Y00). The silicate feat ures as well as the water maser features 
are very stable (Y00; lEngels&Leinertll994l) . so it is natural to as- 
sume that the dust and water are co-spatial, as we can show that 
similar conditions can explain the observations. Both can then be 
locate around the companion in disc, which could have been formed 
while presently C-rich star was still O-rich, and as a result, it con- 
tains the O-rich dust and gas which is responsible for the observed 



O-bearing dust/molecule features. However, a small amount of car- 
bon rich dust from the later stellar wind could also have been cap- 
tured by the companion. We therefore assume that both astronomi- 
cal silicate and amorphous carbon grains are present in the masing 
disc. 

We assume that the dust is optically thin and is heated by the 
car bon star of radius R, = 2 AU and temperat ure of T, = 2400 
K lEneels & Leinertlll994l : benson & Little-MarenirJll996l) . The 
heating of the dust by the companion is negligible if it is a main 
sequence star. The fitting parameters are the dilution factor W, the 
mass of dust, and the dust grain size a. We find that the data are best 
described for a = 0.5 /im. Note, that it is larger than the typically as- 
sume d gra i n size of 0. 1 /an in th e envelopes of AGB stars (see e.g. 
|jud[l996: iGroen ewegen 2006). On the other hand , it seems rea- 
sonab l e that grains would grow whilst in a disc (se£j Ljsjaueretal] 
1 19961 IjuraetaljEoOll : iD'Alessio etafl kOOrl bent et alj l200d) . 
whereas the smallest grains would be blown away by the carbon 
star radiation (see Section l4~TV The observed IR spectrum is well 
described (see Fig. |5J by a model with W = 1.6 10~ 4 for M as = 
5.6 10~ 7 M Q of astronomical silicate and M ac = 6.7 10~ 8 M G of 
amorphous carbon (with corresponding temperatures of T as = 290 
K and r ac = 470 K obtained from the energy balance). The distance 
between the carbon star and the dust material can thus be estimated 
from the dilution factor: D = 0.5r" 2 *, 80 AU. 

3.5 Maser pumping model 

As we argued above, the observed emission of the 615 — 523 maser 
at 22 GHz most probably originates in the disc around the compan- 
ion to the carbon star. The hydrogen and water vapor in the disc 
are probably mixed with the amorphous carbon and astronomical 
silicate dust grains. We assume a solar abundance of elements, giv- 
ing a water-to-ga s mass ratio of /h 2 o = 6 10~ 3 (see, for example, 
IJeong etalJ 2003). The maser strength depends on the hydrogen 
number density n H2 , water fraction / H ,o and the dust-to-gas mass 
ratio f& as well as on the gas and dust temperatures. 

The dust temperature is determined from heating by the car- 
bon star and cooling by its own radiation (see Sect. I3.4i . The gas 
temperature is determined by collisions with the dust grains and 
radiative cooling b y wate r molecules. We take _ 10 8 cm -3 
iKvlafis & Normanl ll99ll) and use / d = 10 2 1 Yamamur a et alJ 
2000 ), a typical ratio in circumstellar envelopes, everywhere in the 
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Figure 5. The silicate emission bands in the spectra of V778 Cyg observed 
by ISO/SWS and IRAS/LRS (Y00). The blue dotted and black dashed curves 
refer to astronomical silicate and amorphous carbon model spectra, respec- 
tively. The dot-dashed curve corresponds to the 2400 K black body emission 
of the carbon star. The red solid curve shows the total model spectrum. 

masing disc. We justify our choice of these parameter values in 
Sect. 14.11 Solving the energy balance for the gas, we find T as 300 K 
at the disc mid-plane. We assume the same temperature throughout 
the disc. 

The conditions in the disc are within the range giving rise 
to collisionally-radiatively pumped, unsaturated maser emission. 
In order to estimate the maser power, we solve simultaneously 
the radiative transfer and statistical balance equations for the first 
45 rotational levels of the ground and (010) excited vibrational 
levels of ortho-H20. We employ the escape probability method 
in order to solve the radiative transfer equation (for details, see 
Babk ovskaia & Poutanenl I2004L 120061) . For a disc half-thickness 
of H ~ 10 14 cm at the outer disc edge (obtained by considering 
the conditions for hydrostatic equilibrium, see below), we obtain 
a maser absorption coefficient of k = -4.5 10~ 14 cm" 1 . This gives 
the maximum maser optical depth r m » 0.73 k R out ~ -20 along 
the line of sight at the disc mid-plane in the direction of the maser 
component at the systemic velocity (see Fig. [5), assuming that am- 
plification takes place only on one side of the disc. The correspond- 
ing maximum brightness temperature is 77, = -T cx exp(-r m ) = 
5 10 lc K, for an excitation temperature T cx = -100 K. Averag- 
ing the brightness temperature over the disc height we get 7b ~ 
1.5 10 10 K. This is above the lower limit of 6 10 8 K imposed by 
the observational data. Thus, the conditions in the disc can easily 
reproduce the observed strength of the maser emission. 



4 DISCUSSION 
4.1 Disc stability 

From the IR spectrum, we have estimated a distance of ~ 80 AU 
between the carbon star and the region containing the dust. This 
separation is consistent with the scenario in which the dust has 



been captured in a disc around a (probably main- sequence) com- 
panion (Y00). Moreover, this estimate is just above the lower limit 
constrained by the velocity stability of the maser spectral features 
I Szczerba et al. 2006). We estimated the radius of the masing disc 
to be about 40 AU by fitting the spatial distribution of masers (Table 
1). This is also consistent with the observed velocity profile. 

For the disc to be stable, its radius should certainly be smaller 
than the distance from the companion to the inner Lagrangian point, 
which is about 44 AU, for the the assumed mass of the carbon star 
of about 1M (Y00) and the companion mass of M s = 1.7M G . 

The stability of the disc also puts constraints on the dust-to- 
gas mass ratio /d. The radiation pressure F ra $ from the carbon star 
acting on the dust should be balanced by the gravitational force 
F grav due to the companion acting on the gas: 

F grav \4nD 2 c j\ R 2 3/ d j < ' 

where Q — 0.2 is the extinction efficiency of the dust (see 
lOssenkopf e7ai]|l992h and p d ^ 3g cm 3 is the dust density. For 
a disc radius of R = 40 AU, equation gives f& < 0.02. 

We use the total dust mass in the disc Ma ugt = M as + M ac = 
5.7 10~ 7 M G to provide an estimate of the hydrogen number density 
assuming a homogeneous isothermal disc: 

«h 2 = f , Mdus ' > 5 10 7 (77300 K)~ 1/2 cnT 3 , (5) 

where Vdi sc — 4HR 2 ut is the disc volume, H = R ou tV t /Vn(R) - 
10' 4 cm is the disc half-thickness (assuming hydrostatic equilib- 
rium), V, m 1.5(77300 K) 1/2 km s _I is the gas thermal velocity, 
V K = 6.1 km s -1 is the Keplerian velocity at the outer disc ra- 
dius R out = 40 AU and mn 2 is the mass of hydrogen molecule. 
An additional condition for disc stability is the requirement that 
the radiation pressure must be smaller than the rate of momentum 
exchange between the grains and the hydrogen molecules, giving 
another constraint on «h 2 : 

F md ^ WH 2 Vt/« co il - mn 2 V^n Hl 7Ta 2 , (6) 

where f co u is the typical time between grain-molecule collisions. 
We obtain wh 2 S; 5 10 7 cnT 3 , which is in perfect agreement with 
Eq. This justifies our choice of n Hl = 10 8 citT 3 and / d = 10~ 2 
when estimating the maser power in Sect. 13. 51 



4.2 Origin of a disc with m = 2 

We note that at the present time there is no theoretical descrip- 
tion of how to produce an m = 2 warped disc (see, how- 
ever, iLubow & Qgilviell2000l and reference therein). An m = 2 
warp might be induced by a non-linear development of warp- 
ing instabilities due to the tidal for ces from the companion 
(Lubow 1992; Lubow & Ogilvie 2000). The m = 1 warping re- 
sults in density waves describ ed by the azimuthal mode m = 2 
jPa p^U2izou & Terqueml H995h . The effect of radiation pressure 
on these waves can produce an m = 2 warp. Such 
m = 2 warps cannot realistically have an amplitude A larger 
than the disc dimensionless half-thickness H/R. Assuming hydro- 
static equilibrium in the disc, we can estimate the ratio H/R = 
Vt/Vjt(R) — 0.24, which is of the same order as the value of 
A x 0.23 which we obtain (see Table 1). More theoretical work 
is required to address the possibilities of m = 2 warps properly. 
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5 CONCLUSIONS 

MERLIN imaging of water masers from the silicate carbon star 
V778 Cyg show clear evidence that they arise from a disc around a 
companion to the carbon star. We estimate the mass of the compan- 
ion to be ~ 1.7 M Q and the distance between the stars to be ~ 80 
AU. This is independently confirmed by the stability over time of 
the velocity of the major maser spectral features. The radius of the 
masing disc is about 40 AU. We develop a model for the disc which 
is consistent with all available observational constraints. 

The silicate dust, observed in the infrared by ISO and IRAS, 
is likely to be co-spatial with the water vapor producing the maser 
emission. This O-rich material must have originated during an ear- 
lier phase of the carbon star and has been captured in a disc around 
the companion. The carbon star is about a thousand times more lu- 
minous than its companion and is the major source of heating in the 
system. 

We find that the observed spectrum of dust radiation is well 
described by a combination of 0.5yum astronomical silicate and 
amorphous carbon grains, with masses M as = 5.6 10~ 7 M o and 
M ac = 6.7 10~ 8 M o . We model the conditions necessary for a stable 
disc, which constrain the dust-to-gas mass ratio to be fa < 0.02 and 
the hydrogen number density to be « H2 <; 5 10 7 cirr 3 . The dust 
temperatures of T w = 290 K and T m = 470 K, obtained from the 
energy balance, yield gas temperature around 300 K (determined 
by collisions with the dust grains and radiative cooling by water 
molecules). These conditions at the same time predict very effec- 
tive water masing and show that the masers could reach a brightness 
temperature more than an order of magnitude larger than the lower 
limit, obtained from the observational data. 

The observed maser geometry and velocity distribution are 
consistent with a doubly-warped Keplerian disc, described by the 
azimuthal mode m = 2. We find that the thickness of the disc 
and the warp amplitude of our model are also consistent with the 
demands of hydrostatic equilibrium. We predict successfully the 
asymmetric appearance of the disc and explain the origin of the 
various maser components. At present, the theory of such m = 2 
discs is absent, and this could be an interesting challenge for future 
theoretical studies. 
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